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Abstract — A new Multi-Carrier Differential Chaos Shift Key- 
ing (MC-DCSK) modulation is presented in this paper. The 
system endeavors to provide a good trade-off between robustness, 
energy efficiency and high data rate, while still being simple 
compared to conventional multi-carrier spread spectrum systems. 
This system can be seen as a parallel extension of the DCSK 
modulation where one chaotic reference sequence is transmitted 
over a predefined subcarrier frequency. Multiple modulated data 
streams are transmitted over the remaining subcarriers. This 
transmitter structure increases the spectral efficiency of the 
conventional DCSK system and uses less energy. The receiver 
design makes this system easy to implement where no radio 
frequency (RF) delay circuit is needed to demodulate received 
data. Various system design parameters are discussed throughout 
the paper, including the number of subcarriers, the spreading 
factor, and the transmitted energy. Once the design is explained, 
the bit error rate performance of the MC-DCSK system is 
computed and compared to the conventional DCSK system under 
an additive white Gaussian noise (AWGN) and Rayleigh channels. 
Simulation results confirm the advantages of this new hybrid 
design. 

Index Terms — Chaos based communication system, Non- 
coherent receiver, Multi-carrier DCSK, Energy efficiency, Per- 
formance analysis. 

I. Introduction 

AS computing devices become ubiquitous, a plurality 
of challenges emerge from the various communications 
paradigms. Some researchers have envisioned that there will be 
"Seven Trillion Wireless Devices Serving Seven Billion People 
by 2020" [1]. In this perspective, spectral and power efficiency, 
interference resistance, security and channel agnosticism are, 
and will continue, to be top requirements for wireless com- 
munication systems. 

Mobile wireless communications performances are dete- 
riorated by device hardware and the propagation environ- 
ments |2|. Fading channels, for instance in Vehicle-to- Vehicle 
communication where moving scatterers cause a non-wide- 
sense stationary uncorrelated scattering (WSSUS) behav- 
ior 0, are typically harsh environments for mobile commu- 
nications. In order to get optimal communication systems in 
varying channels, many techniques can be employed. One is 
the use of multi-carrier systems, such as OFDM, that have 
a high resilience to selective channels if the bandwidth of 
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each subcarrier is smaller than the coherence bandwidth; non- 
coherent communication systems make up the other. It has 
been stated that non-coherent systems can outperform coherent 
ones in fast frequency and time-varying channels, mainly 
because of the spectral inefficiency inherent to the insertion 
of pilots 0]. 

Several combinations of multi-carrier transmission and 
Code Division Multiple Access (CDMA), like Multi-Carrier 
CDMA (MC-CDMA), Multi-Carrier Direct-Sequence CDMA 
(MC-DS-CDMA) and Orthogonal Frequency Code Division 
Multiplexing (OFCDM) are proposed in the literature Q, 
ED, Q. In MC-CDMA, one-bit chips are spread over M 
subcarriers in the frequency domain |5|, while for MC-DS- 
CDMA, time and frequency spreading is used (i.e. TF-domain 
spreading) [7]. Time -domain spreading is employed to increase 
the processing gain in each subcarrier signal, while frequency 
domain spreading is used to increase the total processing gain. 

The chaotic signal has a sensitive dependence upon initial 
conditions property that allows the generation of a theoretical 
infinite number of uncorrelated signals. Those signals have 
been shown to be well suited for spread-spectrum modulation 
because of their inherent wideband characteristic |8| |9| |10|. 
Chaotic modulations thus have similar advantages as other 
spread-spectrum modulations, exempli gratia, including the 
mitigation of fading channels and jamming resistance. The low 
probability of intercept (LPI) performance chaotic signals ifTTI 
agrees with military scenarios and in densely populated envi- 
ronments ( T2). In addition, chaos-based sequences give good 
results as compared to Gold and independent and identically 
distributed sequences for reducing the peak-to-average power 
ratio (PAPR) |[T3l 

A proposed system with a non-coherent receiver, named 
a differential chaos shift keying (DCSK) system, in which 
chaotic synchronization is not used on the receiver side, 
delivers a good performance in multipath channels lfl4ll 
IT31 . Furthermore, differential non-coherent systems are better 
suited than coherent ones for time and frequency selective 
channels [4|. In the DCSK system, each bit duration is divided 
into two equal slots. In the first slot, a reference chaotic signal 
is sent. Depending on the bit being sent, the reference signal is 
either repeated or multiplied by the factor — 1 and transmitted 
in the second slot. A significant drawback of DCSK is that 
half the bit duration is spent sending non-information-bearing 
reference samples 0. This can be accounted as being energy- 
inefficient and a serious data rate reducer. The analytical 
performance derivation of DCSK communication system is 
studied in (16| for fading channels and in IfTTI |[T8l lfl9l lEOl 
for cooperative schemes. The transmission security of DCSK 
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system is improved in ED . In E2l . the spectral efficiency of 
the DCSK is improved, but the system receiver requires an RF 
delay line, which is not easy to implement because of the wide 
bandwidth involved. In a study to overcome the problem of 
RF delay in DCSK systems, Xu et al. proposed a Code Shifted 
Differential Chaos Shift Keying (CS-DCSK) system (23J. In 
their system, the reference and the information bearing signals 
are separated by Walsh code sequences, and then transmitted 
in the same time slot. For such systems, there is no need for 
a delay line at the receiver end. An improved version of the 
high spectral efficiency DCSK system by 11231 is presented 
in (24), where chaotic codes are used instead of Walsh codes, 
with different receiver structures. Because of these advantages, 
some ultrawideband systems based on DCSK or FM-DCSK 
modulations have been proposed for wireless personal area 
networks |25l. E6l. l27l. 

In this paper, we first introduce a new design of a multi- 
carrier DCSK system (MC-DCSK). The system is a hybrid 
of multi-carrier and DCSK modulations. On the transmitter 
side, one of the M subcarriers is assigned to transmitting 
the reference slot, while the other frequencies will carry the 
data slots. In this case, just one chaotic reference is used to 
transmit M — 1 bits, which saves the transmitted bit energy 
and the data rate. After the subcarriers are removed, a parallel 
demodulation is achieved to quickly recover the transmitted 
bits. The proposed system solves the RF delay line problem 
mentioned in [23 1, provides from the properties of DCSK 
system in terms of resistance to interference, increases the 
data rate, and optimizes the transmitted energy of the DCSK 
system with a simple transmitter/receiver design compared to 
conventional multi-carrier spread spectrum systems (i.e MC- 
DS-CDMA). 

Secondly, we thoroughly analyze the performance under 
AWGN and Rayleigh fading channel, without neglecting the 
dynamic properties of chaotic sequences, meaning that in 
our computation approach, the transmitted bit energy is not 
considered as constant. Many approaches have been consid- 
ered for computing the bit error rate performances of the 
DCSK system, such as Gaussian approximation (GA). This 
approximation assumes that the correlator output follows the 
normal distribution. Applied to the DCSK system over an 
AWGN or multipath channel in [28 1, this method provides 
rather good estimates of the BER for very large spreading 
factors, but when the spreading factor is small, the results 
produced by the Gaussian approximation method are rather 
disappointing. Another accurate computation methodology is 
developed in [29], [30] and [31] to compute the BER perfor- 
mance DCSK over different wireless channels. Their approach 
enables the dynamics properties of the chaotic sequence by 
integrating the BER expression for a given chaotic map over 
all possible chaotic sequences for a given spreading factor. 
This latter method is compared to the BER computation under 
the Gaussian assumption in PP . and seems more realistic to 
match the exact BER. However, as indicated in PP . [30|, 
the drawback of the proposed method is the high calculation 
difficulty. Since approaches that have been previously pre- 
sented so far are either invalid for small spreading factors or 
involving highly complex computations, we extend a simple 



and accurate method in this paper for computing the exact 
performance for a single-user MC-DCSK system. The system 
is evaluated first over an AWGN to highlight the problem 
of non constant bit energy. Then, the performance of MC- 
DCSK system is evaluated under Rayleigh fading channel. The 
proposed method includes the computation of the probability 
density function of the chaotic bit energy and the integration of 
the BER over all possible values of the PDF. The advantage 
of this method lies in the fact that it gives an exact BER 
expression without neglecting the dynamical properties of 
chaotic sequences with low computing charges. 

Thirdly, we derive the analytical bit error rate expressions, 
and we show the accuracy of our analysis by matching the 
exact numerical performance, especially when the spreading 
factor is low. We can conclude that the proposed system can be 
suited for Wireless Sensor Network (WSN) applications [32 1, 
which are power-limited and evolve in harsh environments. 

The remainder of this paper is organized as follows. In 
section [Il| the characteristics of chaos-based systems are 
described with an emphasis on DCSK. The III section covers 



the architecture of the MC-DCSK system. The energy and 



spectral efficiencies of the system are examined in section IV 



The performance analysis is explained in section [V] Simula- 
tion results and discussions are presented in section [VI] and 
concluding remarks are presented in section |VII| 



II. DCSK COMMUNICATION SYSTEM 
AND WEAKNESS POINTS 

In this section, the DCSK communication system, together 
with its weakness points, are discussed. The conventional 
DCSK modulation scheme will be used in section |Vl] as 
a comparative to illustrate the performance enhancements 
obtained from the main contribution of this paper. 

A. DCSK communication system 

As shown in Figure [T] within the modulator, each bit Sj = 
{—1, +1} is represented by two sets of chaotic signal samples, 
with the first set representing the reference, and the second 
carrying data. If +1 is transmitted, the data-bearing sequence 
is equal to the reference sequence, and if —1 is transmitted, an 
inverted version of the reference sequence is used as the data- 
bearing sequence. Let 2(3 be the spreading factor in DCSK 
system, defined as the number of chaotic samples sent for 
each bit, where (3 is an integer. During the i th bit duration, 
the output of the transmitter e^fc is 



for 1< k < P, 
for p < k < 2P, 



(1) 



where Xk is the chaotic sequence used as reference and Xk-p 
is the delayed version of the reference sequence. 

Figure [T] illustrates that the received signal r& is correlated 
to a delayed version of the received signal r^+p and summed 
over a bit duration 7], (where 7], = PT C and T c is the chip 
time) to demodulate the transmitted bits. The received bits are 
estimated by computing the sign of the output of the correlator 
(i.e., see Figure [T] (c) the DCSK receiver). 
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Fig. 1. Block diagram of the general structure of the DCSK communication 
system, (a) is the DCSK transmitter, (b) represents the DCSK frame (c) is the 
DCSK receiver. 



B. Weakness of DCSK 

In this paper, we are not working on improving the inherent 
lack of security of non-coherent systems, as the security 
issue was addressed in our previous work [21], where a 
secure chaos-based multi-carrier communication system was 
proposed. Our focus in this work is on the spectral and energy 
efficiencies having a good performance. 

As shown in Figure [T] half the bit duration time is spent 
sending a non-information-bearing reference. Therefore, the 
data rate of this architecture is seriously reduced compared to 
other systems using the same bandwidth, leading to a loss of 
energy. The reference sequence dissipates half the energy of 
each bit. 

III. MULTI- CARRIER DCSK SYSTEM ARCHITECTURE 

The system's architecture is intended to be of low complex- 
ity. Numerous extensions could be performed to this system 
for different performance optimizations. The system presented 
here is in its most elementary form. 

A. Chaotic generator 

In this paper, a second-order Chebyshev polynomial func- 
tion (CPF) is employed 



Xk+i = 1 - 2x fc 



(2) 



This map is chosen for the ease with which it generates 
chaotic sequences and the good performance ll33l . In addition, 
chaotic sequences are normalized such that their mean values 
are all zero and their mean squared values are unity, i.e., 

E(x k ) = and E(x 2 k ) = 1. 

B. The transmitter 

The MC-DCSK system benefits from the non-coherent ad- 
vantages of DCSK and the spectral efficiency of multi-carrier 
modulation. In this system, the input information sequence 
is first converted into U parallel data sequences s u (t) for 



u = 1, 2, ...U. The independent data sequence s u (t) with 
equal probability value is +1, or —1, where 



M-l 



s «(*) = X! s «>*(*) 



(3) 



where s Ut i is the i th bit of the u th sequence data and M — 1 
is the number of data per u th sequence. 

As shown in Figure [2] a reference chaotic code x u (i) to be 
used as a reference signal and spreading code. After a serial- 
to-parallel conversion, the M — 1 bits stream of the u th data 
sequence are spread due to multiplication in time with the 
same chaotic spreading code x u (t). 



K u{t) = ^ x u.kh{t - kT c ), 



(4) 



k=l 



where, is the spreading factor, h(t) is the square-root-raised- 
cosine filter. This filter is band-limited and is normalized to 
have unit energy. Let H(f) = F{h(t)}, where F denotes 
a Fourier transform. It is assumed that H(f) is limited to 
[—B c /2, B c /2] which satisfies the Nyquist criterion with a 
rolloff factor a (0 < a < 1) where B c = (1 + a)/T c . 

The chaotic signal x u (t) modulates the first subcarrier 
as reference, after which the data signals spread by M — 1 
modulate the M — 1 subcarriers. 

Therefore, the transmitted signal of the MC-DCSK is given 
by: 



e(t) = x u (t) cos(27r/it + (j>\)+ 



M 

E 

i=2 



S u i(t)x u (t) COS(27T fit + 4>i) 



(5) 



where <fii represents the phase angle introduced in the carrier 
modulation process. In this paper, we normalize the transmit- 
ted energy in every subcarrier. 

For the MC-DCSK, the modulated subcarriers are orthog- 
onal over the chip duration. Hence, the baseband frequency 
corresponding to the i th subcarrier is f = f p + i/T c , where f p 
is the fundamental subcarrier frequency. The minimum spacing 
between two adjacent subcarriers equals A = (1 + a)/T c , 
which is a widely used assumption Q. 



Reference 



Data 



"*j A=(l+a)/T c I*" 



~*|B c =(l+a)/T c [*~ 



Fig. 3. The power spectral density of a band-limited MC-DCSK system. 

Figure [3] shows the power spectral density (PSD) of the MC- 
DCSK system. Let B be the total bandwidth of the proposed 
system. When both bit duration and B are set, the chip 
duration T c as well as the spreading factor j3 depend on the 
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Fig. 2. Block diagram of the MC-DCSK system 



number of subcarrier M, the bandwidth B c of each subchannel 
or the subcarrier spacing A. In our design, we divide the 
total band B into M equi-width frequency bands, as shown in 
Figure [3] where all bands are disjoint. The bandwidth of each 
subcarrier band B c is: 

B c = (l + a)/T c - 
The total required bandwidth B is: 

B = MB C , 

B = M(l + a)/T c - 
Thus, the spreading factor function of the system parameters 



is: 



P = T b /T c 



T b B 
M(l + a) 



(6) 



In order to validate our system design and evaluate its 
parameters, such as the spreading factor, the bit energy, the 
bit error rate performance, and the spectral efficiency. The 
received signal is: 



r(t) = X(t)e(t) + n(t), 



(7) 



where X(t) is the channel coefficient, n(t) is an wideband 
AWGN with zero mean and power spectral density of Nq/2. 

In our paper, X(t) — 1 for AWGN case. Otherwise, for 
fading channel analysis, has a Rayleigh probability density 
function given by : 



fx(z) = 



z>0- 



(8) 



In our paper, the channel coefficients are assumed to be 
constant during a frame of MC-DCSK symbols and change to 
new independent values from one frame to another. 



C. The receiver 

The block diagram of the MC-DCSK receiver is illustrated 
in Figure [2] One of the objectives of this design was to provide 
a robust receiver providing good performance. We consider a 
set of matched filters, each demodulating the desired signal of 
the corresponding carrier frequency fi, and then the signals 
are sampled every kT c time. The outputs discrete signals 
are stored in matrix memory. The matrix implementation 
of the receiver simplifies the parallel data recovery, where 
the decoding algorithm is described as follow under AWGN 
channel: 

First, at the same time, the output of the first match is stored 
in matrix P and the M — 1 data signals are stored in the 
second matrix S, where: 



P = [x u ,i + n u ,i, x„, 2 + n Ut 2, 



where n Ui k is the k th sample of additive Gaussian noise added 
to the reference signal. 
The matrix S is: 



S= 



s u ,ix u .i + n u 1 



Su.l x u,P + n u,/3 



M-l 



s u ,m-iXu,i3 + n u p 



where n l u k is the k th sample of additive Gaussian noise added 
to the i th bit of u th data sequence. 

Finally, after j3 clock cycles, all the samples are stored, and 
the decoding step is activated. The transmitted M — 1 bits are 
recovered in parallel by computing the sign of the resultant 
vector of the matrix product: 



s u = sign(P x S') 



(9) 



where x is the matrix product and ' is the matrix transpose 
operator. In fact, this matrix product can be seen as a set of a 
parallel correlator where the reference signal multiplies each 
data slot, and the result is summed over the duration j3Tc. 
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IV. Energy efficiency 

The energy efficiency of the proposed system is improved as 
compared to the DCSK system. In fact, for the DCSK system, 
a new chaotic reference is generated for every transmitted bit, 
and in our case, one reference is shared with M — 1 modulated 
bits. For a conventional DCSK system, the transmitted bit 
energy E% is: 



Eh — E, 



data 



+ E re f 



(10) 



where Edata and E re f are the energies to transmit the data 
and reference respectively. Without loss of generality, the data 
and the reference energies are equal: 



E, 



data 



E 



ref 



T 2 
X i,k 



(ID 



k=l 



Then for DCSK system, the transmitted energy Eh for a 
given bit i is: 



E b = 2T C ^ 



T 2 



(12) 



fc=i 



In the MC-DCSK system, one reference energy E re f is 
shared with M — 1 transmitted bit, then the energy of one 
given bit is the sum of its data earner energy and a part of 
the reference energy: 



VP — TP J_ -^ref 
— ^data T - 



(13) 



In our system, the energies on the M subcarriers are equal: 



Edata — E re f — T c ^ ; 



i,k 



k=l 



The bit energy expression function of Edata is: 
M 



E h 



-E, 



data 



(14) 



(15) 



M - 1 * 

To study the energy efficiency, we compute the transmitted 
Data-energy-to-Bit-energy Ratio (DBR): 

Edata 



DBR = 



E h 



For the MC-DCSK system the DBR is: 
DBR M ~ l 



M 



(16) 



(17) 



In a conventional DCSK system (i.e for M — 2), half the 
energy Eh is transmitted into the reference for each bit, and 
then the DBR is: 



DBR 



(18) 



As shown in Figure [4] for M = 2 where we have one new 
reference for every bit bit, in this case, the MC-DCSK system 
is equivalent to a DCSK system with DBR = ~ . This means 
that 50% of the bit energy Eh is used to transmit the reference 
used for one bit. For the same bit energy Eh, in MC-DCSK 




20 3D 4D 50 60 70 

Subcarriers M 

Fig. 4. DBR for a system for various amount of data subcarriers. 



system, we can see for example that for M > 20, the reference 
energy accounts for less than 5% of the total bit energy Eh for 
each bit of the M — 1 data stream. This mean that the energy 
used to transmits the reference is shared with M — 1 bits. 

V. Performance analysis of MC-DCSK 

In this section, the performance of the MC-DCSK system 
is evaluated, and the analytical BER expression is derived 
under an AWGN and Rayleigh fading channels. To derive 
the analytical BER expression, the mean and the variance of 
the observation signal D u i must be evaluated. We start by 
mentioning some properties of chaotic signals which will be 
used later in this paper to compute the mean and the variance 
of the observation. The chaotic generator is very sensitive 
to initial conditions, and we can deduce that the different 
chaotic sequences generated from different initial conditions 
are independent from each other. In addition, the independence 
between the chaotic sequence and the Gaussian noise is also 
true [8]. For the normalized chaotic map with zero mean, the 
variance (Var(.)) is equal to one (Var(x) = E(x 2 ) = 1). 

Since the channel coefficients are assumed to be constant 
and the same for all frequencies during a frame of MC-DCSK 
symbols and change to new independent values from one 
frame to another. The decision variable for the i th bit of the 
u th data stream at the output of the correlator is: 



), (19) 

fe=i 



where X u is the channel coefficient for the u th data stream. 
n u : k and n l u k are two independent zero Gaussian noises 
coming from the reference and the i th bit subcarrier. 



D u ,i — T c S Ut i ^u x u,k + S n u,k n u,k 

k=l k=l (20) 



~\~T C J, X v I A M n^ k + \ u s u,i n u,k J , 
k=l " K ' 
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where T c is the time chip. The i th bit of the u th data stream 
is decoded by comparing the output D u i to a threshold zero. 
For mathematical simplification we set the time chip equal to 
one (Tc = 1). 

In the decision variable equation, the first term is the useful 
signal, while the second and third are zero-mean additive 
noises. 

The output of the correlator for the MC-DCSK of equation by: 
(|20b can be written in the form 



D„ 



(M-i)A^ E, 

: M 



(a) 



BER = 

i Pr ( D u>i < 0| s Utl = +1) + i Pr (D a>i > 0| s Uji = -1) 
= i erfc ( SM^a ^ 

2 C JC \ ^Varl £>„>„<=+!] J ' 

(25) 



where erfc(x) is the complementary error function defined 



erfc(x) 



dfi 



W + Z, (21) The BER expression for the MC-DCSK system is: 



E, 



(«) _ 



M 



M- 1" 



-T, 



E* 2 

k=l 



u,k i 



where E^ is the transmitted bit energy for a given data 
sequence u. 



fc=i 



fe=i 



For a given i bit of an u data stream, the instantaneous 
mean and variance of the decision variable are derived as 
follows : 



BER =-erfc 
2 



MN 



M 2 /3N$ 



(M-l)A^ 2 {M-l)\iE i 



(«)2 



(26) 

Many approaches have been considered for computing the 
BER of chaos-based communication systems, with the most 
widely used being the Gaussian approximation, which consid- 
ers the transmitted bit energy E b as constant H3411 . Because of 
the non-periodic nature of chaotic signals, the transmitted bit 
energy after spreading by chaotic sequences definitely varies 
from one bit to another [35]. Based on this fact, the overall 
BER expression of the MC-DCSK system is then given by: 



BER = 





MN 



(M-1)\ 2 E 2(M-1)A 4 B2 



M &Nq 



f(E b )f(\)dE b d\, 
(27) 



E(D Uji ) 



(M-l)A^ 



(22) 



Since the three terms of pT| ) are uncorrected, the noise 
samples and channel coefficients are independent, the condi- 
tional variance of the decision variable for a given bit i th is: 



Var(D u>i )=E((M SI lxlEl u h u , i r) 



k=l 



E I ( ^u^u.k^u^i^Ujk^) 



k=l 



+E 



(En u ,kK tk ry(^X 2 u E^s u ^ 



Finally after simplification, 

, s (M - l)E { h u) , , , 

Var(Di) - ^ + /3iV 2 /4 ■ 



(23) 



(24) 



In order to compute the BER with our approach, the error 
probability must first be evaluated for a given received energy 
E^ and the channel coeficient X u . Considering the bit energy 
(or chaotic chips) as a deterministic variable, the decision 
variable at the output of the correlator is necessarily a random 



Gaussian variable. Using equations (22i, ( |24| i, the bit error 
probability is : 



where / (A), and p (E^) are the probability density functions 
of channels coefficient and the bit energy E b respectively. 
The BER expression can be simplified as: 




M 



M 2 (3 




(M-l) 7b 2(M-1) 2 7 6 2 



where 7^ = X 2 Eb/N is product of two random variables: 
The bit energy Eb and the channel coefficient A. 

A. BER computation methodology under AWGN channel 

In this section, the AWGN channel is considered to highlight 
the non constant bit energy problem. In this case, the channel 
coefficient is constant equal to one A = 1 and j b = Eb/No. 



In order to compute ( 28 1, it is necessary to get the bit energy 
distribution for a given chaotic map and a spreading factor. 

For this aim, we fitted the histogram of the energy distribu- 
tion for the CPF sequence. Figure [5] shows the histogram of 
the bit energy after spreading by the CPF chaotic sequence for 
j3 = 20. This histogram has been obtained using ten million 
samples. From these samples, energies of successive bits are 
calculated for a given spreading factor. The bit energy is 
assumed to be the output of a stationary random process [36|; 
hence the histogram obtained in Figure [5] can be considered 
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Fig. 5. Histogram of the bit energy distribution for j3 = 20 



as a good estimation of the probability density function of the 
bit energy. 

Given the shape of this energy distribution, the analytical 
expression appears difficult to compute, leaving numerical 
integration as a solution for performing the BER computation. 



The expression ( 28 1 can be computed numerically, taking into 
account the bit-energy variation. 

BER computation methodology under Rayleigh channel 

The objective is to compute the BER given in equation (28 1 
by taking into account both the bit energy and the channel 
coefficient. The analytical derivation of 7^ seems difficult 
because the bit energy have an irregular PDF shape. An 
alternative solution consists to estimate the PDF of 7;,, then 



integrated equation ( 28 i numerically. 



VI. Simulation results and discussions 

To evaluate the effect of the number of subcarriers on 
performance, we plot the computed BER expressions with 
simulation results of the MC-DCSK system over an AWGN 
and Rayleigh channels. The results obtained are for different 
numbers of subcarriers M and spreading factors j3. 

The parameters of the simulation in a mono-user case are set 
as follows: the MC-DCSK system uses the square-root-raised- 
cosine chip waveform for a roll-off factor a equal to 0.25. As 
shown in equation iRS, the spreading factor is computed as 
a function of the number of subcarriers M, the bit duration 
Tf,, and the total allocated bandwidth B. In our simulations, 
we set the bit duration T b = 400, B = 1 and for M = 64 
the allowed spreading factor (3 = 5, for M — 16 subcarriers 
P = 20, for M = 8 P = 40, and f or M = 2 P = 160. 

Figure [6] presents the performances obtained from the BER 
expression of (28 1 and the Monte Carlo simulations of the 
MC-DCSK system under AWGN channel. It clearly appears 
that there are an excellent match between simulations and our 
computed BER expression for any number of subcarrier and 
spreading factor. 

In Figure [7] we study the effect of the number of subcarriers 
on the system performance under AWGN channel. To that end, 




5 10 

Eb/NO [OB] 



Fig. 6. Simulation and BER expression, for different values of spreading 
factor (3, number of subcarrier M, and under AWGN channel 



we set the spreading factor to P = 5 and the bit duration Tb, 
and then we assume that the bandwidth B is wide enough 
to support any number of subcarriers M. Figure [7J shows 
interesting results of our proposed MC-DCSK system in terms 
of performance enhancement. In fact, for a given spreading 
factor, when the number of subcarriers M increases, the DBR 
ratio tends toward one, meaning that less referecne energy 
used to transmit one bit or the reference energy is shared with 
more M — 1 bit. This performance improvement proven in 



expression (28 1 means that for high number of subcarriers M, 
we need less energy to reach a given BER. In this figure, we 
show the performance improvement by simulation for M = 2 
and M = 64, with a fixed spreading factor equal to P = 5. In 
the case of M — 2, the MC-DCSK system is equivalent to a 
DCSK system. The result shown in Figure [7] can be seen as 
a performance comparison between the proposed system with 
that of the conventional DCSK. 

In the same Figure, we can observe a degradation in 
performance between the MC-DCSK system for M = 64 
and the coherent BPSK one. This degradation comes from the 
two noise sources added to the reference and the data carrier 
signals. 

To understand the performance behavior of MC-DCSK 
system for different spreading factor the optimal spreading 
factor must be discussed. The optimal spreading factor was 
studied in ll37l for single carrier DCSK system. Figure [8] 
evaluates the effect of the value of the spreading factor on 
the performance of the MC-DCSK under AWGN channel. 
The simulated bit error rate is plotted for different values 
of spreading factor p with a fixed E^/Nq and number of 
subcarriers M = 2. The bandwidth is assumed to be wide 
enough to support any spreading factor value. Because the 
BER expression is approximated and computed by numerical 
integration, the theory in this case can only qualitatively 
describe the dependence of the MC-DCSK of spreading factor. 
Simulation shows that the intermediate spreading factor values 
between 5 and 50 minimize bit error rates at fixed E^/Nq. 
From this result, we see that good performances are obtained 
for low spreading factor values which makes this system 
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Fig. 7. BER comparison of MC-DCSK for M=64 and DCSK (i.e MC-DCSK 
for M=2) where the spreading factor /3 = 5 under AWGN channel 




—6— Computed BER MC-DCSK, M=2, p = 5, E(J. )=0.7 
— ~ Simulation MC-DCSK, M=2, p = 5, E(l 2 )=0.7 
-^"Computed BER MC-DCSK, M=64, p = 5, E(J. 3 )=0.7 
Simulation MC DCSK, M=6~4, p = 5, E(J. 2 )=0.7 
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Fig. 9. Simulation and BER expression, for a spreading factor j3 = 5, 
number of subcarrier M = 2 M = 64 under Rayleigh channel with average 
gain power E(\ 2 ) = 0.7 



implementation feasible even for moderate bandwidth. 
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Fig. 8. BER versus the spreading sequence length for MC-DCSK for M=2 



Figures [9] and 10 evaluate the effect of the Rayleigh 
channel on the performance of the MC-DCSK system. The 
computed bit error rate given in equation ( |28| and the com- 
puter simulation of the BER are performed. The system's 
performance plotted in Figure [9] is evaluated for two different 
subcarriers M = 2 and M = 64, spreading factor equal to 
j3 — 5, and an average gain power of the channel E(\ 2 ) = 0.7. 



In Figure 10 the performance is evaluated for two different 
subcarriers M = 2 and M = 64, spreading factor equal to 
j3 = 50, and an average gain power of the channel E(\ 2 ) = 
0.9. 

It clearly appears that there are an excellent match between 
simulations and our computed BER expression for any number 
of subcarrier, average gain power, and spreading factor. 

Discussions 

The proposed system meets the following properties: 

• Non-coherent system: Robust receiver; 

• Spread spectrum system: resistance to interferences; 




Computed BER MC-DCSK, M=2, p = 50, E(l )=0.9 
Simulation MC-DCSK, M 2. p = 50, E(J. 2 )=0.9 
Computed BER MC-DCSK, M=64 p = 50, E(i 2 )=0.9 
-Simulation MC DCSK, M=64, p = 50, E(J. 3 )=0.9 



Eb/NO [dB] 



Fig. 10. Simulation and BER expression, for a spreading factor @ = 50, 
number of subcarrier M = 2 M = 64 under Rayleigh channel with average 
gain power E(\ 2 ) = 0.9 



• Chaotic signals: easy to generate, low PAPR in multi- 
carrier transmissions, good correlation properties; 

• Multi-carrier DCSK: high spectral efficiency, low power 
consumption. 

VII. Conclusion 

An energy-efficient non-coherent multi-carrier spread spec- 
trum system has been presented. From the outstanding energy 
inefficiency drawback imposed by time-multiplexed differen- 
tial modulations, a novel frequency multiplexed architecture 
is designed. The multi-carrier characteristic of this novel 
design enables significant energy savings and a higher spectral 
efficiency as compared to differential systems because in the 
new system, the reference signal is only sent once for many 
parallel bits. The energy efficiency of the proposed system is 
analyzed and a DBR is derived, with results showing that 
for M > 20 subcarriers, the energy lost in transmitting the 
reference is less than 5% of the total bit energy per bit. The 
performance of the proposed system is studied, the bit error 
rate expression is derived for Rayleigh fading channels, and 
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simulation results match the theoretical expression, proving the 
accuracy of our approach. To compare the performance of the 
proposed system with that of the DCSK, the simulated BERs 
are plotted with the same spreading factor, where results prove 
an increase in performance as compared to the conventional 
DCSK. Our future work will focus on defining multi-user 
access strategies and performance improvement of this system. 
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